Charge Exchange and Chemical Reactions with Trapped Th 3 



L. R. ChurchillQ M. V. DePalatis, and M. S. Chapman 
School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332-0430 

(Dated: January 4, 2011) 

We have measured the reaction rates of trapped, buffer gas cooled Th + and various gases and 
have analyzed the reaction products using trapped ion mass spectrometry techniques. Ion trap 
lifetimes are usually limited by reactions with background molecules, and the high electron affinity 
of multiply charged ions such as Th 3+ make them more prone to loss. Our results show that reactions 
of Th 3+ with carbon dioxide, methane, and oxygen all occur near the classical Langevin rate, while 
reaction rates with argon, hydrogen, and nitrogen are orders of magnitude lower. Reactions of 
Th 3+ with oxygen and methane proceed primarily via charge exchange, while simultaneous charge 
exchange and chemical reaction occurs between Th 3+ and carbon dioxide. Loss rates of Th 3+ in 
helium are consistent with reaction with impurities in the gas. Reaction rates of Th 3+ with nitrogen 
and argon depend on the internal electronic configuration of the Th 3+ . 

PACS numbers: 34.70.+e,82.30.Fi 



I. INTRODUCTION 

Unlike most atomic nuclei, which have excitation en- 
ergies in the range of kcV to MeV, the nucleus of the 
thorium isotope, 229 Th, has an excited state just several 
eV above the nuclear ground state The transition 
between the nuclear ground state and this unique iso- 
meric state lies within the UV optical spectrum, where 
it can be addressed using coherent light sources. The 
coherent control of the electronic states of atoms with 
tunable lasers has been a major focal point of modern 
atomic physics. Extending this paradigm to the control 
of nuclear states of atoms would represent a significant 
achievement. The transition between the 229 Th nuclear 
states could potentially be used as a frequency reference 
with a fractional uncertainty approaching 10~ 20 [2]. Fur- 
thermore, the transition may be especially sensitive to 
changes in the value of the fine structure constant, al- 
lowing up to 5-6 orders of magnitude enhancement in 
measurements of its time variation 3j. However, this 
latter point requires further study 

In order to demonstrate coherent control of the nuclear 
state, thorium atoms must be confined in such a way that 
they can be continuously interrogated and observed. Un- 
like the lower ionization states, triply ionized thorium has 
a convenient level structure for fluorescence detection and 
laser cooling [2j . The low- lying electronic states of Th 3+ 
and the optical transitions between them are shown in 
Fig. [T] Using the more common thorium isotope, 232 Th, 
our group has demonstrated the creation, trapping, and 
laser cooling of Th 3+ @. 

In our earliest observations of 232 Th 3+ , the ions would 
remain in the trap for only a few seconds. In general, 
elastic collisions with background molecules are not likely 
to result in ejection of ions due to the large depth of an 
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ion trap. The primary modes of Th loss are charge 
exchange and chemical reactions with background gases. 
By improving the background vacuum and increasing the 
purity of the buffer gas used for initial cooling, we were 
able to extend the trap lifetime of Th 3+ to minutes. By 
removing the buffer gas immediately after the initial trap 
loading and laser cooling the ions, a lifetime of > 10 
minutes was obtained ||. 

The limited Th 3+ trap lifetime presents a significant 
challenge in time-intensive experiments like measuring 
the hyperfine states of 229 Th 3+ and searching for the nu- 
clear isomer transition. Furthermore, the short lifetime 
represents a tremendous cost in performing experiments 
with 229 Th 3+ given the extraordinary price of the isotope 
(> $100k/mg). 

To better understand the charge exchange and chem- 
ical reaction processes and quantify their rates, we con- 
ducted a series of experiments to determine the reaction 
rate coefficients between Th 3+ and various gases. A sum- 
mary of the experimental results can be found in Table [I] 
The remainder of this paper is devoted to describing these 
experiments in greater detail. 

A schematic of the ion trap used in these experiments 
is shown in Fig. [T] The trap is outfitted with a chan- 
nel electron multiplier (CEM) for electronic detection of 
trapped ions. When used in combination with the well 
understood mass selection characteristics of a linear ion 
trap, the CEM can be used to determine if a specific 
species of ion is present in the trap. 

The ion trap is loaded with Th 3+ via laser ablation 
of a thorium metal target with the third harmonic of 
a pulsed Nd:YAG laser (A = 355 nm). The target is lo- 
cated near the trap axis and oriented perpendicular to it. 
The voltage on the dc endcap nearest the ablation tar- 
get is gated with the ablation pulse, dropping to ground 
when the pulse is fired, and increasing to 100 V for ion 
confinement. 

Helium buffer gas is present in the system throughout 
the loading process. The ablated Th 3+ ions are initially 
moving at speeds of 10 km/s as they approach the ion 



2 



Reactant 


k 


k s 


(e) 


k L x 10" 9 


He 


< 3.4 x 10" 15 


< 6.0 x 1(T 16 


(0.18) 


1.6 


Ne 


< 4.1 x 1CT 15 


< 1.2 x icr 15 


(0.3) 


1.0 


Ar<*> 


1.3 x 10" 14 


1.7 x 10" 14 


(1.3) 


1.5 


N« 


1.6 x 10~ 13 


1.6 x 10" 13 


(1.0) 


1.8 


H 2 


4.6 x 10~ 13 


2.1 x 10~ 13 


(0.46) 


4.4 


CH 4 


2.6 x 10" 9 


3.6 x 1(T 9 


(1.4) 


2.8 


o 2 


3.7 x 10~ 9 


3.7 x 1(T 9 


(1.0) 


1.7 


C0 2 


2.8 x 10" 9 


4.0 x 10~ 9 


(1.4) 


1.8 



TABLE I: Charge exchange and chemical reaction rate coef- 
ficients for Th 3+ . Reaction rate coefficients are given in units 
of cm 3 s _1 . Here, kL is the calculated Langevin rate, and k 
is the experimental reaction rate coefficient. The corrected 
reaction rate coefficient is k £ — e ■ k, where e is the ion gauge 
correction factor for the gas in question. Experimental reac- 
tion rate coefficients for helium and neon are upper bounds. 
The accuracy of the coefficients given here is limited by the 
accuracy of the ion gauge used to measure the pressure of the 
reactants. (*) The reaction rate for Ar and N2 depends on 
the electronic configuration of Th 3+ . The experimental rate 
refers to reaction from the 5F 5 / 2 ground state. 
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FIG. 1: (Color online) (a) Th 3+ energy levels and optical 
transitions, (b) Schematic drawing of ion trap system. A 
target outside of the ion trap (left) is ablated using the third 
harmonic of a pulsed Nd:YAG laser (diagonal line). Lasers 
used for optical excitation of the trapped ions are aligned 
along the trap axis (horizontal line). A CEM located outside 
of the trap (right) is used for electronic detection. The CEM 
is mounted off-axis to allow clear optical access. 



trap. As the hot Th 3+ ions collide with room tempera- 
ture helium atoms, they lose some of their thermal en- 
ergy. This process serves two critical purposes. First, 
it reduces the energy of some ablated Th 3+ ions as they 
traverse the length of the ion trap, thereby increasing the 
likelihood they will be trapped. Secondly, it damps the 
motion of trapped Th 3+ ions, allowing them to come to 
thermal equilibrium at some fraction of the trap depth. 
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FIG. 2: (Color online) Accumulated fluorescence signal for 
determination of Th 3+ loss rate. This data was taken with 
5 x 10 -5 torr of helium in the vacuum chamber. The loss rate 
was determined to be r' 1 = 0.0100(2) s _1 . 



The helium pressure is commonly held between 10~ 5 and 
10~ 4 torr in these experiments. Generally, the number 
of ions loaded increases with higher buffer gas pressure, 
while the final equilibrium temperature of the trapped 
ions decreases. 



II. HELIUM 

The loss of Th 3+ from the ion trap in 5 x 10~ 5 torr 
of helium buffer gas is shown in Fig. [5J Here, the Th + 
decay rate is measured by observing the reduction in flu- 
orescence over a series of images. This method can be 
used only when the lifetime of the ions is sufficient to 
allow enough fluorescence measurements for a proper ex- 
ponential fit. 

When the Th 3+ trap lifetime is too short for observa- 
tion of fluorescence, mass selective electronic detection 
with the CEM is employed to measure the decay rate. 
We load ions into the trap via ablation and wait for some 
amount of time. The rf and dc voltages applied to the 
trap are then tuned such that only Th 3+ remains stable. 
The contents of the trap are then delivered to the CEM 
by lowering the voltage of the trap endcap nearest it. 
The time interval before the mass selection is varied, and 
the data is accumulated over several iterations for each 
time interval. Many iterations are required to suppress 
the variation due to fluctuations in the number of ions 
loaded. This method is typically used when the Th + 
lifetime is a few seconds or less. 

The decay rate through any given reaction channel is 
proportional to the density of the reactant in the system. 
This relation can be written as r _1 = kn, where n is the 
density of the reactant, and k — (av) is the reaction rate 
coefficient for the given channel. Here, a is the reaction 
cross-section and v is the relative velocity between the 
reactants. For a given reaction, the rate coefficient can 
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FIG. 3: (Color online) Th 3+ loss rate as a function of BIP 
helium buffer gas pressure. The slope of the linear fit implies 
a reaction rate coefficient of k = 3.4 x 10 -15 cm 3 s _1 (fc E = 
6.0 x 10~ 16 cm 3 s _1 with an ion gauge correction factor of e — 
0.18), while the intercept implies a background loss rate of 
t" 1 = 3.3 x 10" 3 s" 1 . It is likely that the Th 3+ reacts with 
an impurity in the helium rather than with the helium itself; 
therefore, the value for k sets an upper bound on the reaction 
rate coefficient for helium. 



function of the pressure of the reactive gas. The reaction 
rate coefficient is extracted from the slope of a linear fit 
of this data. The intercept of the linear fit gives the loss 
rate due to background gases in the vacuum. For helium, 
the data and corresponding linear fit are shown in Fig. [3] 

Helium is a noble gas, and a charge exchange reaction 
between it and Th 3+ is endothermic by 4.6 eV. This sug- 
gests that such a reaction is strongly inhibited for low 
energy reactants. Therefore, the loss shown in Fig. [2] is 
most likely due to an impurity in the helium rather than 
to the helium itself. 

The helium gas we use in our experiment comes from a 
cylinder with a built-in purifier (BIP). The BIP helium is 
supplied by Airgas, Inc. According to the specifications, 
the BIP helium gas has an impurity level < 1 ppm. The 
data of Fig. [3] implies a reaction rate coefficient between 
the BIP helium and Th 3 + on the order of 1CT 16 to lfr 15 
cm 3 s _1 . Thus, the reaction between Th 3+ and the con- 
taminant within the BIP helium responsible for its loss 
proceeds with a reaction rate coefficient on the order of 
1(T 10 to 1CT 9 cimV 1 . 

A reaction rate coefficient of 10~ 9 cm 3 s _1 is consistent 
with the value predicted by a classical collision model 
introduced by Langevin @ and refined by Gioumousis 
and Stevenson [lipj . This model is helpful in estimating 
the reaction cross-sections between ions and molecules. 
We assume the ion is a point particle with charge Ze, 
and that the molecule is spherically symmetric with po- 
larizability a. The primary interaction between the ion 
and molecule is due to the electric field of the ion and 
the field-induced diploe moment of the molecule. The 
interaction potential scales as r~ 4 . 

Below threshold values of the relative velocity and the 
impact parameter, the interaction potential leads to a 



spiraling orbit of the ion and molecule about their shared 
center of mass. If a reaction between the ion and molecule 
is exothermic, it is assumed to occur with unit probability 
when a spiraling orbit occurs. The so-called Langevin 
cross-section for reaction is then [Toj 
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where v is the relative velocity between the particles and 
/x is the reduced mass. The corresponding reaction rate 
coefficient is 



kt = 1ixZe\foni 



(2) 



The preceding equations are given in Gaussian (cgs) 
units, which is how they are most commonly found in 
the literature. 

The Langevin model is not a comprehensive theoretical 
formulation, and therefore, it only suffices to provide an 
estimate of the reaction rate between a given ion and re- 
actant molecule. The model does not apply to molecules 
that possess a permanent dipole moment. In that case, 
the leading term in the interaction potential scales as 
r~ 2 rather than r~ 4 [ll[. Furthermore, several experi- 
ments [l2l-[l4| have shown a departure from the inverse 
velocity dependence of the classical cross-section in some 
reactions. Ultimately, the probability of a reaction occur- 
ring in the event of a classical interaction is dependent 
on the combined energy state of the reactants, the avail- 
able energy states in the final products, and the possible 
reaction pathways between them 15|. Several experi- 
ments [16h20| |. including our own, have found that some 
reaction rates can differ significantly depending on the 
internal energy state of a reactant. 

Nevertheless, the Langevin rate is useful in setting an 
upper bound on the reaction rate between ions and non- 
polar molecules [2l|. Many ion- molecule reactions do 
occur at or near the Langevin rate [l7], I22U26I ] . The 
Langevin rates for the reactions we investigated are given 
in Table U with the experimental results. 



III. GASES WITH LOW REACTION RATES 

Loss rates of Th 3+ in the presence of neon and ar- 
gon are shown in Fig. [U while loss rates in the presence 
of nitrogen and hydrogen are shown in Fig. [5] For all 
of these decay measurements, a partial pressure of 10 -5 
torr of helium was added to the system to cool the ions 
for fluorescence imaging. As can be seen from Fig. [3J 
this quantity of helium does not add significantly to the 
background loss rate of Th 3+ . The reactions with argon, 
nitrogen, and hydrogen all proceed relatively slowly, with 
reaction rate coefficient orders of magnitude below the 
corresponding Langevin rates. Since our interest was in 
gases that react with Th 3+ near the Langevin rate, we 
did not attempt to identify the products of these reac- 
tions. 
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FIG. 4: (Color online) Th 3+ loss rate in the presence of 
neon and argon. A partial pressure of 10" 5 torr of he- 
lium was added to the test gas to cool the ions for fluores- 
cence observation. The reaction rate coefficient for neon is 
k = 4.1 x KT 1S cm 3 s _1 (fe e = 1.2 x 10~ 15 cm 3 s _1 with an 
ion gauge correction factor of e = 0.30). The decay rates 
shown for argon represent Th 3+ loss from the 5F 5 / 2 ground 
state. The reaction rate coefficient is k — 1.3 x 10~ 14 cm 3 s -1 
(fc e = 1.7 x 10 -14 cm 3 s _1 with an ion gauge correction factor 
of e= 1.3). 



to rf micromotion could conceivably contribute to over- 
come an otherwise endothermic energy gap |27i-l29j , such 
as with helium and neon. However, such micromotion- 
enhanced reactions would require micromotion ampli- 
tudes much larger than is likely. While some contri- 
bution of the micromotion to the loss rates is possible 
for endothermic reactions, our results are consistent with 
quoted impurities [30j | . 

Argon is also a noble gas, but a charge exchange reac- 
tion between it and Th 3+ is exothermic by 4.2 eV. The 
Th 3+ reaction rates with argon and nitrogen were found 
to depend on the electronic configuration of the Th 3+ . 
The data in Fig. [4] represent loss rates from the 5F 5 / 2 
ground state of Th 3+ . The 5F 5 / 2 state was isolated for 
study by shuttering the 690 nm laser between fluores- 
cence measurements. This allows Th 3+ in the A-manifold 
(see Fig. [T|) to be pumped back into the ground state by 
the 984 nm laser. By keeping the time between measure- 
ment sufficiently long, the systematic effect introduced 
during measurement was held below 10%. 

IV. Th 3+ EXCITATION EFFECTS ON 
REACTION RATES 
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FIG. 5: (Color online) Th 3+ loss rate in the presence of hy- 
drogen and nitrogen. A partial pressure of 10 -5 torr of he- 
lium was added to the test gas to cool the ions for fluores- 
cence observation. The reaction rate coefficient for hydrogen 
is k = 4.6xl0~ 13 em's"' (fc e = 2.1 x 10" 13 cm 3 s _1 with an ion 
gauge correction factor of e = 0.46). The decay rates shown 
for nitrogen represent Th 3+ loss from the 5F 5/ / 2 ground state. 
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Like helium, neon is a noble gas with a high ioniza- 
tion energy that inhibits charge exchange reactions with 
Th 3+ . Such reactions are in this case endothermic by 1.6 
eV. As can be seen in Fig. 21 Th 3+ loss rates in neon were 
never significantly larger than what could be attributed 
to the presence of the helium buffer gas. In addition, the 
reaction rate coefficient we determined from the data is 
consistent with reaction occurring between the Th 3+ and 
0.5 ppm impurity in the gas. 

It is worth noting that the kinetic energy of Th 3+ due 



State dependent effects on electron capture in colli- 
sions between neutral atoms and multiply c harg ed ions 
have been studied extensively (see, e.g., [3ll434l |). The 
effect of Th 3+ optical excitation on its reaction rate with 
argon can be seen clearly in Fig. [S] In one set of mea- 
surements, the 690 nm laser was shuttered between fluo- 
rescence measurements. This is the same data set shown 
in Fig. SI In the other two sets of measurements, the 690 
nm laser was left on continuously. Increasing the power 
of the 984 nm laser increases the population of the 6D 5 / 2 
excited state, and hence the reaction rate, until the ex- 
cited state is saturated. The corresponding saturation of 
the argon and nitrogen reaction rates are shown in Fig. [7J 

If the 5F7 /2-6D5 / 2 manifold is approximated as a two- 
level system, the excited state population obeys the re- 
lation 
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where lsat,eff is the effective saturation intensity. The ra- 
tio of the effective saturation intensity to the natural sat- 
uration intensity, / sat = 2-7r 2 /icr/3A 3 , is approximately 
equal to the ratio of the Doppler-broadened linewidth, 
Yd, to the natural linewidth, Y. Therefore, 
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(4) 



By fitting the data in Fig. [7] to a function in the form 
of Eq. ([3]), we find an effective saturation intensity of 
isat.eff ~ 40 mW/cm 2 . According to Eq. (gj), this corre- 
sponds to a Doppler-broadened linewidth of Yd ~ 300 
MHz, which is consistent with an independent measure- 
ment of the transition linewidth. 
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FIG. 6: (Color online) Effect of atomic excitation on Th 3+ 
loss rate in the presence of argon. The data that was taken 
without 690 nm light represents loss from the Th 3+ 5F 5 ; 2 
ground state. Increasing the power of the 984 nm laser in- 
creases the reaction rate by increasing the number of ions in 
the 6D5/2 excited state, from which reaction is more likely. 



The reaction rate from the 5F 7 / 2 electronic state was 
found to be equal to within experimental error to the 
reaction rate from the 5F 5 / 2 state for both argon and ni- 
trogen. The 5F7/2 state was isolated for measurement by 
shuttering the 984 nm laser between measurements, while 
continuously applying the 690 nm laser. In this way, the 
long-lived 5F7/2 state was continuosly repopulated. To 
further demonstrate that optical excitation to the 6D 5 / 2 
state was responsible for the increased reaction rates, we 
verified that detuning the lasers far from resonance had 
the same effect as shuttering them. 

By optically exciting with the 1087 nm laser, we found 
that reactions with nitrogen were also faster from the 
6D3/2 state than from the 5F 5 / 2 state. No quantitative 
comparisons were made between the reaction rates from 
the 6D3/2 state and the 6D 5 / 2 state. The reaction rate 
between Th 3+ and argon with 1087 nm optical excitation 
was not measured. 



V. GASES WITH HIGH REACTION RATES 

Loss rates of Th 3+ in the presence of carbon dioxide, 
methane, and oxygen are shown in Fig. [5] All of these 
gases react relatively quickly with Th 3+ at rates com- 
parable to the Langevin rate. Even at partial pressures 
below 10 -8 torr, reactions occurring in the space between 
the ablation target and the trapping region significantly 
reduced the number of Th 3+ ions loaded. To ensure suffi- 
cient loading and an adequate SNR for fluorescence mea- 
surement, 10~ 4 torr of helium was added to the system 
for each measurement. The loss rates measured here were 
in all cases significantly higher than the loss rate with 
only 10~ 4 torr of helium present. 

The products of reaction between Th 3+ and carbon 
dioxide, methane, and oxygen were identified using mass 
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FIG. 7: (Color online) Saturation of Th 3+ reaction rate in 
the presence of (top) 3 x 10~ 6 torr argon and (bottom) 3 x 
10~ 6 torr nitrogen. In each case, 1 x 10~ 5 torr helium is also 
present. 



selective CEM detection. For these measurements, 2 x 
10~ 8 torr of the test gas and 10~ 4 torr of helium buffer 
gas was introduced into the system. The ion trap was 
loaded a number of times. Each time the trap was loaded, 
we waited a short period (~ 1 s) for reactions to occur, 
performed a mass selection, and checked for the presence 
of ions. Since the rf power supply for our trap could not 
provide high enough voltage to perform rigorous mass 
selection on singly ionized molecules, we were able to 
identify only doubly ionized molecules. 

When either carbon dioxide or oxygen was present in 
the system, Th 2+ and Th0 2+ were found in the trap. 
Although thorium dioxide is chemically stable, no Th0 2 + 
was detected with either gas. In a previous work with 
singly-ionized thorium, Johnsen et al. [22| found that 
oxidation occurs via the sequential reactions 
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They found that the first reaction proceeded quickly, with 
a rate coefficient of 6 x 10 -10 cm 3 s _1 , while the rate of 
the second reaction was more than an order of magnitude 
lower, with a coefficient of 2 x 10~ n cm 3 s _1 . 

Reactions between Th 3+ and methane, CH4, resulted 
in Th 2+ and ThCH 2 + . Special care was taken to prop- 
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FIG. 8: (Color online) Th 3+ loss rate in the presence of (a) 
carbon dioxide, (b) methane, and (c) oxygen. Measurements 
were made with 10 -4 torr of helium buffer gas. The data point 
at zero pressure is the loss rate when only helium at 10 -4 torr 
is present. Reaction rate coefficients for each case are listed 
in Table U 

erly identify ThCHj + . Andrews and Cho 2(y\ were able 
to create ThCH4 by ablating thorium in a methane en- 
vironment. However, Margalo et al. 35] found that the 
only reaction channel between singly-ionized thorium and 
methane is 

Th+ + CH 4 ->■ ThCHj + H 2 . (7) 

Only 1 amu/e separates the mass-to-charge ratio of 
ThCH 2+ from ThCH 2+ and Th0 2+ . By carefully map- 
ping the voltage space in which Th0 2+ and ThCH 2 + 



FIG. 9: (Color online) Doubly-ionized thorium molecules. 
Ions are identified by comparing the voltage space in which 
they are stable to the theoretical stability region, (a) The 
reaction between Th 3+ and carbon dioxide produces Th0 2+ . 
(b) The reaction between Th 3+ and methane produces Th 2+ , 
which in turn reacts with methane to produce ThCH?, + . (c) 
The reaction between Th 3+ and oxygen produces Th 2+ , which 
in turn reacts with oxygen to produce Th0 2+ . 



were stable, we were able to achieve the resolution nec- 
essary to correctly determine the product (see Fig. [5]) . 

Once the reaction products were identified, we were 
able to measure how Th 3+ evolved over time in the pres- 
ence of each of these gases. The results are shown in 
Fig. [TU] for carbon dioxide, methane, and oxygen. Each 
data point in these graphs represents the average and 
standard deviation of six measurements. For each mea- 
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FIG. 10: (Color online) 232 Th 3+ in the presence of (a) carbon 
dioxide, (b) methane, and (c) oxygen. The data is scaled 
according to the values given in Table [TTJ The data sets are 
fit to the system of equations given in Eqs. (fT51) - ([T5)) . 



surement, the trap was loaded, and Th 3+ was mass se- 
lected. A fluorescence measurement would be taken im- 
mediately after the mass selection. After the specified 
wait time, the ion of interest would be mass selected and 
the contents of the trap delivered to the CEM. The result- 
ing CEM signal was normalized according to the initial 
fluorescence measurement. By normalizing the signal in 
this way, the impact of loading fluctuations on the data 
was minimized. 

As was mentioned previously, the rf power supply for 
our trap could not provide high enough voltage to per- 
form rigorous mass selection on singly ionized molecules. 



However, by ramping the rf to its maximum voltage, we 
could isolate ions with mass-to-charge ratios greater than 
156 amu/e. The only ions in the trap above that mass-to- 
charge ratio would be singly-ionized thorium molecules. 
These are represented as ThX + in Fig. [TO] 

Direct comparison of the CEM signals of the various 
ions is of limited value since the gain of a CEM can vary 
from ion to ion. Ions with higher charge states experience 
a greater acceleration in the field created by the high 
negative bias at the entrance of the dynode. These ions 
strike the dynode surface with greater energy, resulting in 
a higher initial emission of electrons. The CEM gain also 
depends on mass, as lighter ions tend to cause a larger 
response than more massive ones. The chemical nature 
of the ion can also play a role. 

We verified several times that the Th 3+ CEM signal 
was linear with measurements of total fluorescence over 
the dynamic range of both instruments. However, the 
specific slope of the relationship at any given time was 
highly sensitive to the background light level around the 
system. Curtains were employed to stabilize and mini- 
mize background light; however, small shifts in the cur- 
tains incurred while the system was prepared for data 
collection still had a noticeable effect. 

In order to compare the relative populations of the 
various ions in the trap over time, we adopted a simple 
method for scaling the CEM data. The method is based 
on two assumptions. We assume first that ions are not 
lost from the trap over the period of investigation. This 
is reasonable given that the trap depth is on the order of 
100 eV for the ions in question, while the exothermicity of 
the reactions is only on the order of a few eV. We further 
assume that the only ions present in the trap are those 
identified in our earlier search for the reaction products. 
Since our method of searching for reaction products fo- 
cused only on products resulting from fast reactions, this 
assumption could lead to error when considering longer 
time scales. However, here we focus our attention on the 
rapid chemical kinetics that occur on short time scales. 

To scale the data we consider the sum 



N(t) = ^a irH (t), 



(8) 



where N is the total trap population, and ni and an are 
the fluorescence normalized CEM signal and the CEM 
scaling factor for ion i, respectively. The sum is over all 
ion species present in the trap. The scaling factors are 
determined via a least-squares algorithm that attempts 
to make N(t) = 1 for all t. The data shown in Fig. QJJ] 
was scaled using this method. The scaling factors are 
given in Table QTJ 

It is clear from Fig. [TUJ that the dynamics responsi- 
ble for the appearance of Th0 2+ in the presence of car- 
bon dioxide are different from the dynamics that cause 
Th0 2+ to arise in oxygen and that bring about ThCH?, + 
in the presence of methane. Based on the identified prod- 
ucts, the possible reactions between Th 3+ and these gases 
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Species 


Carbon Dioxide 


Methane 


Oxygen 


Th 3+ 


9.0 


4.6 


3.6 


Th 2+ 


23.4 


26.2 


14.7 


Th0 2+ 


36.8 




12.7 


ThCH 2 ^ 




53.9 




ThX+ 


20.5 


30.0 


39.3 



TABLE II: CEM scaling factors for TV+ reaction data 



can be written generally as 



Th 3 



Th 



3+ 



Th 3 



co 2 



CH 4 



2 



02, 



71 

72, 



2+ 



Th 
ThO 2 



CO^ 



CO H 



Th 2+ + CH+ 



ThCHo 



2+ 



Th : 
ThO 







2 + 



H 2 

2 

0+ 



(9) 
(10) 

(11) 



Here, a,, /3j, and ji represent the relative probabilities of 
the reaction channels. The possible subsequent reactions 
with Th 2+ can be summarized as 



Th 



2 + 



Th 



2 + 



co 2 



CH 4 



Th 



2 + 



o 2 





ThX+ 4 






Th0 2+ - 




03 


ThX+ + 




<3 4> 


ThCH| + 


+ H 2 


73 


ThX+ + 




74^ 


Th0 2+ H 


-0 



(12) 
(13) 
(14) 



No decay in the quantity of Th0 2+ or ThCH 2 ^ is seen 
over the time scales investigated, so these can be treated 
as stable final products. Thus, from the above reaction 
equations, we can write a system of differential equations 
describing the chemical kinetics in each gas. For example, 
in carbon dioxide, 



( Th3+ ) 
dt v ' 

|„(ThO-) 
-|„(Th*+) 



-km (Th 3+ ) (15) 
-k 2 n (Th 2+ ) + aihn (Th 3+ ) (16) 
a 2 fcin (Th 3+ ) + a 4 k 2 n (Th 2+ )(17) 
a 3 k 2 n (Th 2+ ) (18) 



Similar systems of equations can be written to describe 
the reaction dynamics in methane and oxygen. 

To estimate the relative probabilities of the various 
reaction channels, the data of Fig.fTOlwas numerically fit 
to dynamic equations of the form shown above. Fits were 
performed in the order of the above equations. Once a 
value was obtained for a fitting parameter, that value was 



Branching Ratio 


Carbon Dioxide 


Methane 


Oxygen 




0.36 


0.89 


1.0 




0.78 


0.0 


0.09 




0.10 


0.0 


0.08 


[a, P,l\ A 


0.44 


0.94 


0.83 



TABLE III: Branching ratios for thorium reactions. 



Gas 


Minimum Purity 


co 2 


o 2 


Hydrocarbons 


Argon 


99.9997% 


< 0.5 


< 0.2 


< 0.2 


Hydrogen 


99.999% 


< 0.5 


< 1 


< 0.5 


Neon 


99.999% 


< 0.5 


< 0.5 


< 0.5 



TABLE IV: Purity specifications for research grade argon, 
ultrahigh purity hydrogen, and research grade neon used in 
reaction experiments. Values listed are in ppm. All gases 
were purchased from Airgas, Inc. and specifications are listed 
in the Airgas 2009 Product Catalog [30t | . 



enforced on subsequent fits. The branching ratios that 
were determined in this fashion are given in Table Mil 

The values given in Table IIIII can only be considered 
rough estimates of the actual branching ratios because of 
the uncertainties inherent in the CEM scaling procedure. 
They are given only to illustrate the general reaction 
dynamics. While reactions of Th 3+ with methane and 
oxygen proceed predominantly via charge exchange, si- 
multaneous charge exchange and chemical reaction is the 
dominant branch in reactions between Th 3+ and carbon 
dioxide. Since methane has the highest ionization energy 
and lowest polarizability among hydrocarbons (36j, the 
charge exchange reaction rates between other hydrocar- 
bons and Th 3+ are likely as high as those of methane. 

It is worth considering whether impurities of these fast 
reacting gases could explain the Th 3+ losses in the pres- 
ence of argon, nitrogen, and hydrogen. C0 2 , 2 , and 
CH4 each react strongly with Th 3+ at rates very close 
to the calculated Langevin rates. In order for these 
molecules to account for the observed loss rates in the 
other gases, they would have to be present as impurities 
at the 10~ 4 level reacting at the Langevin rate, consider- 
ably higher than the levels specified by the gas supplier 
(see Table [Vj) . Furthermore, none of these molecules ex- 
hibited different reaction rates due to optical excitation 
of Th 3+ , so their presence as impurities could not explain 
the reaction rates observed for argon and nitrogen. For 
these reasons, it is unlikely that impurities of the type 
measured to react at the Langevin rate are responsible 
for the losses observed for the slowly reacting gases. 



VI. CONCLUSION 

We have determined the effect of various gases on the 
trap lifetime of Th 3+ . Reactions involving carbon diox- 
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ide, oxygen, and methane proceed at a rate near the clas- 
sical Langevin limit. Although charge exchange reactions 
with them are exothermic, the reaction rate coefficients 
for nitrogen, hydrogen, and argon are orders of magni- 
tude less than the Langevin rate. The reaction rate co- 
efficient between Th 3+ and helium provides a measure 
of the impurities in the buffer gas. The measurement of 
the coefficient here provides a standard against which the 
vacuum and buffer gas quality of future systems can be 
measured. 
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